function A=rotateAboutX(phi)
  A=[1,0,0;
     0,cos(phi),-sin(phi);
     0,sin(phi),cos(phi)];
